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Abstract 

We postulate the second-order derivative equation with four pa- 
rameters for spin-1/2 fermions possessing two mass states. For some 
choice of parameters fermions propagate with the superluminal speed. 
Thus, the novel tachyonic equation is suggested. The relativistic 20- 
component first-order wave equation is formulated and projection op- 
erators extracting states with definite energy and spin projections are 
obtained. The Lagrangian formulation of the first-order equation is 
presented and the electric current and energy-momentum tensor are 
found. The minimal and non-minimal electromagnetic interactions of 
fermions are considered and Schrodinger's form of the equation and 
the quantum-mechanical Hamiltonian are obtained. The canonical 
quantization of the field in the first-order formalism is performed and 
we find the vacuum expectation of chronological pairing of operators. 

1 Introduction 

The study of models which allow the superluminal behavior of particles is 
of definite theoretical interest. Previously, models of superluminal particles 
(tachyons) were considered in pQ, [2] (see refs. therein). Here, we pay at- 
tention on the possibility of the superluminal propagation of fermions in the 
framework of the higher derivative (HD) theory. HD models attract atten- 
tion due to their improved renormalization properties [3], [I], [5], [6], [7]. We 
formulate here the second-order (in derivatives) equation for spin-1/2 fields, 
which admits the superluminal particle speed and is the generalization of 
Barut's equation [8], [9] describing fermions with two mass states. 

The paper is organized as follows. In Sec. 2, we formulate the second 
derivative equation for free fermions with the dispersion relation allowing 
the superluminal propagation. The first-order relativistic wave equation is 
derived and the relativistically invariant bilinear form is obtained in Sec. 3. 
The mass and spin projection operators are given in Sec. 4. We obtain the 



electric current and energy-momentum tensor in Sec. 5. Sec. 6 is devoted to the 
introduction of minimal and non-minimal electromagnetic interactions. The 
Schrodinger form of the equation is given and quantum-mechanical Hamilto- 
nian is obtained in Sec. 7. We consider the canonical quantization of a model 
and obtain the the vacuum expectation of chronological pairing of operators 
in Sec. 8. In Sec. 9 the novel tachyonic equation is suggested. The results are 
discussed in Sec. 10. Appendices A and B contain some useful products of 
matrices. We use Euclidian metrics and notations as in [10] , and the system 
of units h — c — 1 is chosen. Greek letters run 1,2,3,4, and Latin letters run 
1,2,3. 



2 The Model 

Let us investigate the second-order (in derivatives) field equation with four 
parameters describing spin-1/2 particles. We suggest HD equation 

7A - -j^ + M + b llll ,d ll + M o75 ) V(aO = 0, (1) 

where d u = d/dx v = (d/dx m ,d/d(it)), i/j(x) is a Dirac spinor and the 
Dirac matrices 7 M obey the commutation relations 7^7^ + 7i , 7/J = 25^, 
75 = 7i 7 2 7 3 7 4- The parameters a and b are dimensionless, and M, M 
have the dimension of the mass. The suggested equation violates in general 
the parity conservation. Equation (1) is the generalization of the Barut's 
equation [8], [9J. At b = 0, M = 0, one comes to the Barut equation which 
was interpreted as the equation for a description of e and fi leptons (see also 
[TT]). One can obtain from Eq.(l) the tachyonic first-order Dirac equation 
by putting a = 0, 6 = 0, M = (o/M = 0). We pay also attention here on 
the particular case of Eq.(l) when parameters a, b, M and M are connected 
with three mass parameters Ki, k 2 , k 3 as follows: 

M = — a = , b = — , M = 6«i. (2) 

At this case Eq.(l) is reduced to the equation 

(flvdv + K i) + «a + K375) H x ) = °- ( 3 ) 

At K3 = 0, one arrives at Barut's equation [8], |9j. In general, four parameters 
a, b, M, and Mq in Eq.(l) are independent. In momentum space Eq.(3) 
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becomes 

(ip + Ki) (ip + n 2 + k 3 7 5 ) ij){p) = 0. (4) 

where p = p^^, p^ = (p,ipo). After multiplying Eq.(4) by the matrix 
(ip — k 2 + K375) (ip — Ki), one obtains the dispersion equation 

(f + 4) (p 2 + 4- 4) = 0, (5) 

where p 2 = p 2 — p 2 ,, and po is the energy of the fermion. It follows from Eq.(5) 

that fermions possess two mass states with masses K\ and At 
k 2 > k 3 fermions have real masses, but at k 2 < k 3 the mass k becomes 
imaginary that indicates on the presence of a tachyon. In this case (k 2 < K3) 
fermions propagate with the superluminal speed. If the mass Ki is huge, such 
state of fermions is not observable. One can speculate that such state of fields 
with big mass K\ and weakly interacting with fields contribute to dark matter. 
In the case k 2 = 7^ 0, we have massless state. Another possible way to 
realize a massless state is to choose n 2 = k 3 = 0. If k 3 = i<p, where <p is a real 
number, the second state is not tachyonic and corresponds to the mass k = 
\J k 2 + ip 2 . Such situation is realized in the the Nambu— Jona-Lasinio model 
with CP-violating condensate ip = ifyj^ip) 7^ [12]. Thus, introducing 
three mass parameters allow us to investigate differing mass states: massive, 
massless and tachyonic. Propagator corresponding to the model based on 
Eq.(4) is given by 

[(ip + ki) (ip + k 2 + k 3 7 5 )] _1 
_ (ip -K2 + K375) (ip ~ K i) 

K 2 — K 2 
where Kj — /S^2 ^3" 

Thus, we have the difference of propagators in Eq.(6) 
so that the state with the mass K\ is a ghost. Another possibility is to 
interpret the state with the mass k as a ghost. These two cases correspond 
to different signs in the Lagrangian. We interpret here the state with the mass 
Kj \ clS db ghost state. When the mass K\ — > 00 that state is not observable. 
But ultraviolet behavior of the propagator with the mass k is improved by 
this additional subtraction. The situation is similar to the Pauli— Villars 
regularization. The difference is that the state with the mass K\ is in the 
spectrum. It should be noted that similar scheme with higher derivative 
equations was explored in the Lee— Wick model [6]. In the general case, when 
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four parameters a, 5, M, and M Q are independent, the dispersion equation 
corresponding to Eq.(l) reads 

oV + M 2 (l + 2a - b 2 ) p 2 + M 2 (M 2 - M 2 ) = 0. (7) 

This equation is valid for a particular case M = 0, b ^ neglecting Eqs.(2). 
Then finding the roots of Eq.(7) at Mq = 0, one obtains squared of masses 
of the field states 

1 + 2a - b 2 + yj(l + 4a - b 2 ) (1 - b 2 ) , 

(8) 

1 + 2a - b 2 - ^(1 + 4a - b 2 ) (1 - b 2 ) . 

If 6 = ±1, one has m\ = m\ and for a > fields have subluminal speed but for 
a < fields propagate with superluminal speed and are tachyons. At b = 0, 
a ^ 0, we arrive at the case considered in [llj. If |6| < 1 and 1 + 2a — b 2 > 
fields have subluminal speed but for 1 + 2a — b 2 < fields propagate with 
superluminal speed. Thus, Eq.(l) allows us to consider different models of 
superluminal particles. 



2 _ 1V1 

1 ~ 2a 2 




3 First-Order Field Equation 

Let us introduce 20-component function which is the direct sum of a bispinor 
ip{x) and a vector-bispinor ip^x): 

»(«) = = ( *w j , (9) 

where the index A runs values A = 0, /i. The vector-bispinor is given by 

M*) = -jj d M x )- ( 10 ) 

Thus, the wave function $?(x) presents the reducible representation of the 
Lorentz group. Now, we use the elements of the entire matrix algebra e A ' B 
[13] with matrix elements and the product of two matrices: 

(e M > N ) AB = S MA 6 NjB , e M > A e B > N = 6 AB e M > N , (11) 
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where indexes A, B, M, N run five values 0, 1, 2, 3, 4. Any matrix with ma- 
trix elements Q>ij Cclll be written as a^e 1 '- 7 , where we imply a summation on 
repeated indexes. Then Eq.(l), taking into account Eqs.(9)-(11), becomes 

d v (V>° + ae^ + (1 + b l5 )) AB * B (x) 

r ,„ i (12) 

+ [M + + M e°'°7 5 J AB * B (a;) = 0. 

Equation (12) can be cast into the first-order relativistic wave equation 

(v v d u + m + M r 5 ) = o, (13) 

where 20-component matrices are given as follows: 

r„ = (e"»° + ae ^) ® J 4 + ® 7l/ (1 + 6 7 s) , T 5 = © 7s , (14) 



and I4 is the unit 4x4 matrix acting in the bispinor subspace, © is the direct 
product of matrices. At b = 0, M = Eq.(13) is converted into equation 
considered in [TT]. To investigate the properties of first-order relativistic wave 
equation (13), one can use the general methods described in [Hj. Eq.(13) is 
convenient for different applications as the matrices of the equation, T u , r$, 
are expressed through the elements of the entire matrix algebra (11). 

The generators of the Lorentz group in the representation of the wave 
function (9) 

[(1/2, 0) © (0, 1/2)] © {(1/2, 1/2) © [(1/2, 0) © (0, 1/2)]} 
are given by [11] 



•V = - ^) ®h + h®\ (Wu - 7,7m) > (15) 
and obey the commutation relations: 

[Jfiu, Jap] — &vuJy,fi + fififiJua ~ &vfiJ pa ~ d^aJvfi. (16) 

The relativistic wave equation (13) is form-invariant under the Lorentz trans- 
formations. Indeed, one can verify that matrices (14) obey the necessary 
commutation equations 

[Ia, Jfiu] = - dxu^ti, [r 5 , J^A = 0. (17) 
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The relativistically invariant bilinear form is 

vj7# = y+rjty, (18) 

where is the Hermitian-conjugate wave function and the Hermit ianizing 
matrix rj is given by [H] 

v = (as m > m - ae iA - e°'°) <8> 74, (19) 

?7 + = 77, and obeys equations: 

rjTi = -rtv, ^r 4 = r+r/, ^r 5 = -r 5 ^ (2 = 1,2,3). (20) 

It should be noted that the Lorentz-invariant is imaginary, ('ITs'I') 4 " = 

— (tyT^). Therefore, to formulate the variational principle, one has to imply 
that the parameter M is imaginary, M = iM\, where Mi is a real parameter. 
Then the "conjugated" equation reads 

(r^V M — M — iM 1 T 5 ) = 0. (21) 

Eq.(21) follows from Eq.(13) (at Mq = iM\) by Hermitian conjugation and 
multiplying the equation by rj. 



4 Mass and Spin Projection Operators 

Here, we restrict our consideration by putting M = (Mi = 0); then masses 
of fermionic fields mi and m 2 are given by Eq.(8). Because there are two 
field states, one can introduce an additional quantum number r in such a 
way that r = 1, 2 for states with masses mi and m 2 , correspondingly [15J. 
Let us consider states of particles with definite energy, p Q , and momentum, 
p. For the definite mass state, we have 

Po = \fp 2 + mj, or p = yp 2 + m\. (22) 

The index r in four momentum p = (p, ipo) is omitted here. Solutions to 
Eq.(13) with definite energy-momentum in the form of plane waves are given 
as follows: 

¥£(x) = -±=U(±p) ex V (±tpx), (23) 
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where V is the normalization volume. Substituting Eq.(23) into Eg. (13), we 
obtain (at M = 0) 

(±ip + M) U(±p) = 0, (24) 

where p = T^p^, and the solution ^>^(x) corresponding to the positive en- 
ergy p describes particles and the solution \&L p ^(x) with negative energy — p 
corresponds to antiparticles. Similar to QED [ID] , we use here the normal- 
ization condition 

/ ^(x)T,¥^(x)d 3 x = l, (25) 

where ^(x) = (^g(x)) + r?, and the integration over the volume V is im- 
plied in Eq.(25). Eqs.(23)-(25) are valid for two values of energy-momentum 
corresponding to additional quantum number r = 1,2. The 20-component 
functions U(±p) obey equations: 

U{±p)T,U{±p) = -2ip„ (26) 

2 

U(±p)U(±p) = T 2^, (27) 

where r = 1 corresponds to p 2 = —m 2 , and r = 2 corresponds to p 2 = — m|. 
Eqs.(26),(27) are similar to normalization conditions for bispinors in QED. 
To obtain the projection matrix extracting solutions to Eq.(24), we use the 
minimal equation for the matrix p = p^T^ (see Appendix A): 

j5 5 - (l + 2a - b 2 ) p 2 p 3 + a 2 p 4 p = 0. (28) 



With the help of the method of the work [16J, one obtains from Eq.(28) the 
projection matrices 

±ip (M T ip) [P 2 - (1 + 2a - b 2 ) p 2 - M 2 ] 
± 2M 2 [(1 + 2a - b 2 ) p 2 + 2M 2 ] ' 1 ' 

that obey equations 

(±ip + M) n ± = 0, (30) 

= n±, n+n_ = o. (31) 

The operator of the spin projections on the direction of the momentum p is 
given by [11] 

i 

° v = ~ 91 — ~\ tabc PaJbc-, (32) 
^IPI 



where |p| = y pi + p% + P3, and obeys the matrix equation 

< - 1) {"l ~ i) = »■ < 33 > 

From Eq.(33), we find the spin projection operator [TT] 

^ = ^ U ± I) U - (34) 



2V~2j\ p 4 
which obeys the equation 

<r P P±i/2 = ±2 P±1 /2- ( 35 ) 

One can verify that the operators p and a p commute [p, <7 P ] = and, therefore, 
they have the common eigenfunctions. Using Eqs.(30),(35 ), one obtains the 
projection operators for pure spin states in the form of matrix-dyads 

U s (±p) • U s (±p) = iV±P ±1/2 n ± , (36) 

with matrix elements (u s (±p) • U s (±p)j = (U s (±p)) A (u s (±p))^, and 
N± is the normalization constant. We have introduced the spin index s = 
±1/2 in Eq.(36). Eq.(36) can be used for calculations of different quantum 
processes of particles obeying Eq.(l) (for M = 0) in the first-order formalism. 
We use the normalization for pure spin states as follows: 

2 

U s {±p)U r {±p) = TSsrjj, U s {p)U r {-p) = 0. (37) 

Eq. (27) follows from Eq.(37) after a summation over the spin index s = 
±1/2. Taking into consideration the relationship (see Appendix A) 



P [Pl/2 + P-l/2) = P, (38) 

one can obtain from Eqs.(36),(37) an expression for matrix density for impure 
spin states 

X; U s {±p) • U 3 (±p) = N ± U ± . (39) 

s 

Using the equation trll± = 2 (see Appendix A) and taking the trace in both 
sides of Eq.(39), we find 

Y,U s {±p)U 3 {±p) = 2N ± . (40) 



Comparing Eq.(40) with Eq.(27), one obtains 

Eqs.(29),(39) generalize expressions obtained in [TF] on the case 6^0. 



(41) 



5 Electric Current and Energy- Momentum Ten- 
sor 

Now, we consider the case Mo = iMi, then Eqs.(13),(21) follow from the 
Lagrangian 



C 



^{x) (T u d u + M + iMtTs) V(x) 



(42) 



-#(z) (r M B„ — M — i M 1 T 5 ) *{x) 
which is a real function. The electric current density is given by [18 

dC dC 



d(d^{x)) d(d^(x)) 
From Eqs.(42),(43),(9),(ll), we obtain 

j^x) = M(x)T^(x) = 



(43) 



(44) 
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iipix^f, (1 + 675) ip(x) + — Yl){x)d^il)(x) - [d^{x)j ip(x) 



Expression (44) does not depend on the mass parameter Mi and includes 
as the usual Dirac current (at b = 0) as well as Barut's convective terms. 
One can verify with the help of Eq.(l) (at M = zMi) that the current 
(44) is conserved: d^j^x) = 0. The charge density follows from Eq.(44) 
and is jo = — iji — ^{x)T^{x) which vanishes for neutral particles. The 
canonical energy-momentum tensor, obtained from the standard procedure 
[IB] , is given as follows: 



T, 



(45) 
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We have taken into account here that Lagrangian (42) vanishes for fields 
obeying the field equations (13) (for M = iM x ) and (21). From Eqs.(9),(ll),(45), 
one finds the energy-momentum tensor 

T, u = ^ip(x) lfl (1 + 675) d v iP(x) - 1 (dj${x)) 7m (1 + &7s) 

a r - - - (46) 1 

+ 2M ( d ^^J d ^{ x )-{ d ^ d ^{ x )) tp(x)+(d v il;(x)J d fJr ip(x)-ip(x)d v d IJ ,ip(x) , 

that is conserved: dyT^ = 0. The canonical energy-momentum tensor (46) 
can be symmetrized by the Belinfante procedure. The energy density is 
given by £ = T 44 . At 6 = 0, M\ = 0, one comes to expressions of j M and 
obtained in [TT] . 



6 Non- Minimal Electromagnetic Interactions 

In this section, we consider the general case with four independent parame- 
ters, a, b, M, and Mq. The minimal electromagnetic interaction is introduced 
by the substitution <9 M — > P M = <9 M — ieA^, where is the four- vector po- 
tential of the electromagnetic field. We also consider the non-minimal elec- 
tromagnetic interaction by adding terms with two parameters ko, ki char- 
acterizing fermion anomalous electromagnetic interactions. The first-order 
relativistic wave equation with minimal and non-minimal interactions reads 
(see also |TT] for a case 6 = 0, M = 0) 



T^D„ + - {k P + k x P x ) IV> + M + M r 5 



V(x) = 0, 



(47) 



where P = e°'° ® h, Pi 



I4 are the projection operators so that 



the relations P 2 = P , P 2 = Pi, P + Pi = 1 hold. We have introduced the 
matrix 



[IV 



r r 



r r 



= (l - 6 2 ) ® (7m% - 7,7m) + * (e*" - e"'") (8) J 4 
+ (>° - ae ^) ® 7 „ (1 + 675) - (e"'° - as '") ® 7^ (1 + 675) • 



(48) 



From Eq.(47), one can obtain equations for Dirac's spinor (bispinor) ^ and 
vector-spinor ip^ with the help of Eqs.(9),(ll),(48): 



lu (1 + 6 75 ) D u + ik 1 - b 2 ) j^jvF^ + M + M o7 5 ^(x 
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(49) 

+ [aD^ + ik a^y v (1 + 675) 7 VV \ ip^x) = 0, 

{Dp + iki^vF^v) ip(x) + (M5^ u + ikxaJF^) ip u {x) = 0, (50) 

where = d^Ay — dyA^ is the strength of the electromagnetic field. Solving 
the matrix equation (50) for vector-spinor ip^ and replacing it in Eq.(49), we 
can find an equation for bispinor if). Such equation includes as minimal as well 
as non-minimal electromagnetic interactions of fermions within our model. 
To clear up the physical meaning of the parameters ho, k\, one can consider 
non-relativistic limit of Eqs.(49),(50). We leave this for further investigations. 



7 Schrodinger Form of the Equation 

Let us consider minimal interactions of fermions by putting ko = k\ = 
in Eq.(47). To obtain the Schrodinger form of the equation and quantum- 
mechanical Hamiltonian we use the method of the work [19] . Thus, we rewrite 
Eq.(47) as follows: 



iT A d t m(x) 



T a D a + M + M r 5 + eA r 4 



The matrix obeys the equation as follows (see Appendix B): 



T\ - (1 + 2a - b 2 )Tl + a 2 A = 0, 



where the projection operator A (A 2 = A) is given by 

A = ( £ 0,0 + £ M) 3 h 



(51) 



(52) 



(53) 



The matrix A extracts the 8-dimensional sub-space of the dynamical compo- 
nents of the wave function (9): $(x) = A\I/(x). To separate the dynamical 
and non-dynamical components of the wave function ^(x), we introduce the 
projection operator II (II 2 = II): 



n 



A = e" 



14, 



(54) 



which defines non-dynamical components Q = Il\E'(x). Let us introduce the 
matrix 



(l + 2a-6 2 )r 4 -r 



(55) 
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so that BY A = A. Multiplying Eq.(51) by the matrix B and taking into 
account Eq.(55), one obtains 



1 r 

+ ^7 



id t $(x) = eA $(x) 



;i + 2a - 6 2 )r 4 - I* (T a D a + M + M T 5 ) *(x). 



(56) 



From the equality A + II = 1, we find the relation = + Q(x). 

One can eliminate the non-dynamical components Q(x) from Eq.(56). In- 
deed, multiplying Eq.(51) by the matrix II, with the help of equations (see 
Appendix B) nr 4 = 0, nT 5 = 0, Iir a IT = 0, we obtain 



Iir a D a $(x) + MQ(x) = 0. 



(57) 



Taking into account the equality ^(x) = &(x) + Q(x), and replacing £l(x) 
from Eq.(57) into Eq.(56), one finds the Schrodinger form of the equation 

id t $(x) =n®(x), (58) 

where the quantum-mechanical Hamiltonian is given by 

U = eA + \ [(1 + 2a - b 2 )T 4 - F\ 
a A L 

x (r a D a + m + M r 5 ) (i - ^nr 6J D 6 

The 8-component wave function <&(x) in Eq.(58) presents only dynamical 
components which describe fermionic fields with two mass states. The Hamil- 
tonian (59) can be written in the simple form using the elements of the entire 
matrix algebra (see Appendix B) 



(59) 



U 



* M 
eA + — 
a 



as 



0,4 



I 4 + e 



4,0 



75 + - 

a 



.4,0 



7 4 (1 + 675) 

<8 7m(l + 675)] An - ^ (^ 4 '° ® h) D 

Eq.(58) with the help of Eq.(60) takes the form of two equations 

id t ip(x) = eA ip(x) + Mtp^x), 



a 



(60) 
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M 

— 7 4 (1 + 075) 
a 



i>4 (x) 



(61) 



+ - 



a 



M + M o75 + 7m(l + ^)D m - —D 



ip(x) 



Eqs.(61) can also be obtained from Eqs.(49),(50), at ko = ki = 0, after the ex- 
clusion of non-dynamical (auxiliary) components ijj m {x) = —(l/M)D m i/;(x). 
So, only components with time derivatives enter Eqs.(61) and Eq.(58). In 
some applications the Schrodinger form of an equation has advantage because 
it contains only dynamical components. In addition, the matrix Hamiltonian 
(60) possesses only non-zero 8x8 matrix block because the wave function $ 
has 8 components. 



8 Field Quantization 



In this section, we consider the case Mq = Mi = but b 7^ 0. Then there 
are three independent parameters a, b, M, and Eqs.(8) are valid. Here, we 
follow to the work [T7] very closely. The general solution to Eqs.(13),(21) (at 
Mq = Mi = 0) in the second quantized theory can be written as follows: 



E M£ 



+), 



X) 



(62) 



E 



a T , P , s U T:S (p) exp(ipx) + b+ U TtS (-p) exp(-ipx) 



*r(x)=E 



a+ s ¥rts\x) + b Ti M7(x) 



(-)/ 



E 



1 



(63) 



a+ PjS U T)S (p) exp(-ipx) + b TtPtS U T)S (-p) exp(ipx) 



w V^PoV 

for the values of the quantum number r = 1,2 corresponding to two mass 
states (22). We introduce the creation and annihilation operators of parti- 
cles a+ ps , a T)P)S , and antiparticles b* PtS , b TjPtS . They obey the commutation 
relations 

{P"r,p,si p' s'l = ^ss'^tt'^pp'j { a r,p,si Or',p',s'} = 0, {tt T _-, Qy/ „/ s / } = 0, 

{K,p,s, K',p',s'} = Sss'Kr'Spp', {&r,p, S , K> ,p> ,s>} = 0, {K >p>s , K',p',s'} = > ( 64 ) 
{a T ,p,s,&r'y,s'} = { a T,p,s,K',p',s>} = { a t,p,s^r',p',s'} = { a t,p,s ' K> ,p> ,s> } = - 
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The energy density of particles and antiparticles is given by 

£ = T 44 = ^(x)T 4 d ^(x) - % - (d V(x)) T^(x). (65) 

With the help of the normalization condition (25) and Eqs.(62)-(64), we find 
from Eq.(65) the energy operator 

H = / £d 3 x = Po (at,p,s a T, P ,s ~ W,p, s K,p,s) ■ ( 66 ) 

T,p,S 

One can obtain from Eqs.(62)-(64) the commutation relations of fields for 
different times t, t'\ 

{*tM(«), ^tn(x')} = {^ t m(x),^ tN (x')} = 0, (67) 

{V t m(x),%n(x')} = K tMN (x,x'), (68) 
K t mn{x, x') = K^ MN (x, x') + K tMN (x, x ), 

KmnM = E (^^) N > ( 69 ) 

We find from Eqs.(62),(63) the functions [TO] 



k ?mn( x >^) = E^T/ ( f/ ^(p)) A f (^r, s (p))^exp[ip(x-x')], 



(70) 



With the help of Eqs.(39),(41), and taking into account the relation p 2 = 
—m 2 , one finds from Eqs.(70): 

( ip(M- ip) (p 2 + Cm 2 T - M 2 ) m 2 T \ 
KMN{X) = ? I 4^(2^ -Cm?) J MJV 6XP( ^ } 

( YjJ^ (M - r„3„) [Cm^ - M 2 - (TM 2 } m 2 T \ ^ 1 ^ 
" 2M3 (2M^ - Cm?) \ 2 Po V 6XPlWj ' 

, , ^ f ip (M + ip) (p 2 + Cm 2 T -M 2 )m 2 \ . . . 
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(72) 

' TjJ^ (M - Y v d v ) [Cm 2 — M 2 — (T^) 2 ] m 2 T \ 1 

2M*(2M 2 -Cm 2 ) ) MN p 2p V 6 %pX) > 



where C = 1 + 2a — b . Exploring the singular functions [10] 

A +( x ) = J2w-T7 ex ^( i P x )^ A _(x) = ^TTexp(-zpx), 
p ^Po v p ^Po v 

(73) 

Ao(x)=i(A+(x)-A-(x)), 
from Eqs.(69),(71),(72), we obtain 

K M f ( T A (M - T u d u ) [Cm 2 T — M 2 — (I^f] m 2 \ 

K ™ n{x) = - % [ 2M3 (2M 2 - Cm?) J MJV ° ( } ' 

(74) 

Eqs.(67)-(74) are valid for two indexes r = 1,2 corresponding to two mass 
states (22). One can verify with the help of Eq. (28), that functions K~(x), 
K+(x) obey the equations as follows: 

(Tpd^ + M) K~(x) = 0, (Tpd^ + M) K+(x) = 0. (75) 

The vacuum expectation of chronological pairing of the operators (the prop- 
agator) is defined by [TO] 



K M {x)K N {y) = k c tM n{x - v) 

(76) 

= 6(x - y Q ) K+ MN (x -y)-9(y - x ) K~ MN (x - y), 

with theta-function 6{x). Then using of Eqs.(71),(72), we find the vacuum 
expectation of chronological pairing of the operators for two values r = 1,2: 

Km{x)Kn{v) 

(77) 

A c (a; - y), 



(M - T v d v ) [Cm 2 T — M 2 — (TM 2 } m 2 ' 



2M 3 (2M 2 - Cm 2 ) 



MN 



where 

A c (x -y) = B{x Q - y ) A+(x - y ) + e(y - x ) A_(x - y). (78) 
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In Eq.(77), we have two propagators corresponding to r = 1 and r = 2 
(there is no summation on index r in Eq.(77)). Thus, equation (77) obtained 
generalizes the expression for the vacuum expectation of chronological pairing 
of the operators [17] on the case of b ^ 0. The density matrix (36) and the 
vacuum expectation of chronological pairing of the operators (77) can be used 
for calculations of different quantum-electrodynamics processes of fermions 
within the model considered. The anticommutator (68) for equal times can 
be obtained from Eq.(74) using the properties of singular functions [10J. 



9 The Novel Tachyonic Equation 

Now, we consider a particular case of a tachyonic equation at b = ±1, M = 0, 
a < 0. Thus, Eq.(l) becomes 

7m (1 ± 75) d, - + M] ${x) = 0. (79) 

This equation describes tachyons with the mass squared (8) at b = ±1, 
m 2 = M 2 /a. In this case mass states degenerated and the field possesses 
only one mass state. As a result, the propagator can not be presented in 
the form of the two propagators difference (6), and the problem of the ghost 
presence can be avoided. It was shown in [7] that in HD theory with the 
degenerated mass states there are not problems with unitarity. As the mass 
of the tachyon depends on two parameters M and a, one can choose for 
simplicity a = — 1. Then using b = —1, we arrive at the novel equation 



ip(x) = 0. (80) 



describing tachyons with the mass squared m 2 = —M 2 . Thus, the tachyonic 
mass in Eq.(80) is M. One can choose another value of b = 1 in Eq.(79). 
The first term in Eq.(80) contains the projection operator (1 — 7s)/2 which 
extracts the left (chiral) component of the bispinor ip(x). The Lagrangian 
corresponding to the first-order formulation of Eq.(80) is given by Eq.(42) 
at Mi = with the matrices (14) with b — — 1, a — — 1. As was already 
mentioned, the Dirac first-order tachyonic equation does not admit the La- 
grangian formulation and, therefore, the theory of tachyonic fermions sug- 
gested has an advantage. 
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For superluminal particles with the tachyonic mass M the energy and 
momentum are given by 



M Mv 

Po = 7^f> P = 7^T' (81) 

so that Pq = p 2 — M 2 and v > 1. From Eq.(81), one obtains 




2 



v 2 -l= - . (82) 



Eq.(82) indicates that if the mass of a tachyon is constant the speed of a 
particle increases with decreasing the energy. For energies much larger than 
the tachyonic mass M, the particle speed approaches the speed of light. From 
Eq.(82), we find the relation 



M = Po Vv 2 - 1 ~ Po V26, (83) 

where the deviation from the speed of light is 5 — v — 1 <C 1. Eq.(83) 
allows us to calculate the tachyonic mass M if the deviation 5 is measured. 
The novel theory of tachyon particles based on Eqs.(79),(80) allowing the 
Lagrangian formulation is of definite theoretical interest. 



10 Conclusion 

We have analyzed the model based on HD Eq.(l) with four parameters a, b, 
M, Mo that generalizes the Barut's model [8], [9] with two parameters a and 
M. This allows us to consider different cases with superluminal propagation 
of fermions. Thus, when parameters a, b, M, M are connected by relation 
(2) with /t 3 > k 2 fermions possess the superluminal speed. Another case of 
superluminal propagation of fermions is realizes when Mo = 0, |6| < 1 and 
a < {b 2 — 1)/2. We pay attention here on a special chose of the tachyon equa- 
tion when b = ±1, a = — 1. Thus at b = — 1, a = — 1 the novel equation (80) 
describes tachyons with one mass state, m\ = m 2 = —M 2 . Because mass 
states degenerated there is no conflict in this case with unitarity (see [7]). We 
have demonstrated that the first-order formulation of equations obtained is 
convenient for finding the conserved electric current and energy-momentum 
tensor, Schrodinger's form of the equation and quantization. The density 
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matrix (36) and the vacuum expectation of chronological pairing of the oper- 
ators (77) allow us to calculate different quantum processes of fermions in the 
framework of the model based on HD equations. We also introduced two pa- 
rameters ko and k\ characterizing non-minimal electromagnetic interactions 
of fermions. There are some questions about the interpretation of parameters 
ko, ki introduced that requires further investigations. The Schrodinger form 
of the equation is convenient because it contains only dynamical components 
of the wave function. In the Schrodinger picture the wave function possesses 
only 8 components for describing fermionic field with two mass states in the 
first-order formalism. 

It should be noted also that if M is a real value and fermions propagate 
with the superluminal speed there is a difficulty of the Lagrangian formula- 
tion of the model because the value ^ITs 1 !/ is imaginary. The same situation 
occurs in the tachyonic Dirac equation when a = 0, 6 = 0, M = (a/M = 0) 
in Eq.(l) as the value ipj^ip is imaginary. This indicates that the tachyonic 
first-order Dirac equation is not healthy. But if M = %M\ (Mi is real), no 
difficulties with the Lagrangian formalism. If four parameters are restricted 
by Eq.(2) and Mq = %M\ there is subluminal propagation of fermions. At 
M = 0, we still have the opportunity to consider the superluminal propaga- 
tion of fermions at the case |6| < 1, a < (b 2 — l)/2 or if b = ±1, a < and 
there are no difficulties with Lagrangian formulation. 

Appendix A 

Now, we find expressions for matrices entering the first-order wave equa- 
tion (13). From Eq.(14), one obtains 

p = p v V v = J(°) ® p (1 + b l5 ) + /W (g) J 4 , (Al) 
where p = Pu'ju, and 

/( o) = £ o, 0) = Vv ( e „,o + a£ o,^ ? /( 2) = ViiVvE w. (A2) 

With the help of Eq.(ll), it can be verified the equations as follow: 

j(0)2 = j(0) j(0)j(l) + j(1)j(0) = j(l) j(2)2 = p 2j(2) ) 

jiO)^) = 7 (2) 7 (0) = Q) 7 (2) 7 (1) = p 2 7 (l) 7 (0) ) 7 (1) 7 (2) = p 2 / (0) / (l) j (A3) 
7 (2) 7 (1) + 7 (D 7 (2) = p 2 7 « 7 (1)2 = a ^2 /( 0) + /( 2)^ 
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Using Eqs.(ll),(A.3), we obtain matrices: 

p 2 = (l + a - b 2 ) p 2 I^ <g> I 4 + ® P (1 + &7s) + a/ (2) ® h, (AA) 

p 3 = (l + 2a - b 2 ) p 2 I (0) ®p(l + b l5 ) + (l + a - o 2 ) p 2 I (1) <g> h 
+a/ (2) ® £(1 + 075), (A5) 
j5 4 = p 2 (l + 2a - b 2 ) p 2 -a 2 p 2 (p 2 /(°) ® J 4 + /( 2 ) <g> J 4 ) , (A6) 

p 5 = (l + a - 6 2 ) (l + 3a - 6 2 ) p 4 /( ) <g> p (1 + 675) 

+ (l + 3a + a 2 - 3a6 2 - 26 2 + 6 4 ) p 4 I {1) <g> J 4 

+a (l + 2a - 6 2 ) p 2 / (2) ® p (1 + 675) . (A.7) 
The matrices (A.1)-(A.7) obey the "minimal" matrix equation: 

p 5 - (l + 2a - 6 2 ) p 2 j5 3 + a 2 p 4 p = 0. (A8) 

From Eqs.(A.l)-(A.8), we obtain traces of the matrices: tr(p)=tr(p 3 )=tr(j5 5 ) = 
0, tr(p 2 ) = 4p 2 (l + 2a - 6 2 ), tr(p 4 ) = 4p 4 [(l + 2a - 6 2 ) 2 - 2a 2 ]. Eq.(A.8) can 
be used for different calculations within the first-order formalism. 

Appendix B 

In this Appendix, we obtain matrices which can be used in the Schrodinger 
form of the wave equation. With the help of Eqs.(ll),(14), we find matrices 

F 2 - 
1 4 — 



F 3 

1 4 



(1 + a - 6 2 )e°'° + ae 4 ' 4 ] ® h + (e 4 >° + ae°> 4 ) <g> 7 4 (1 + 075) , (B.l) 
(1 + a _ & 2 )£ o,o + a£ 4,4] 074(1 + 67s) 

+(1 + a - b 2 ) (e 4 >° + a£°' 4 )(g) 7 4 , (B.2) 

T 4 = { [(1 + a - 6 2 ) 2 + a(l - 6 2 )] + a(l + a - b 2 )e 4 ' 4 } ® 7 4 

+ (1 + 2a - 6 2 ) (e 4 - + as ' 4 ) ® 74 (1 + 075) , (5-3) 
From Eqs.(B.l),(B.3), one obtains the equation as follows: 

T 4 - (1 + 2a - b 2 )T\ + a 2 (5°'° + 5 4 ' 4 ) = 0. (BA) 
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In in the Schrodinger equation, we use matrices: 



a 1 



[l + 2a-b 2 )T 4 -T\ 



1 / a n n/i\_ T 1 



BT m = -e ifi ® 7m (1 + b l5 ) + e^ m ® J 4 , 5r 5 = -e 4 ' ® 7s , (£.6) 
a a 

nr m = e mfi <g> i 4 , flr m Ajir n £> n = e 4 - <g> / 4 z^. (s.7) 

Also relations 

5n = o, r 5 n = o, e 4 ' m $ = o (5.8) 

hold. 
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